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Abstract 

This  work  examines  the  diffusion  of  rubidium  through  a  small  diameter  tube  alone 
and  in  the  presence  of  noble  gases.  A  fluid  dynamics  analysis  is  investigated  to  aid 
in  choosing  a  method  for  transferring  atomic  rubidium  vapor  that  is  reliable  and 
efficient.  Solutions  to  the  time  dependant  ordinary  differential  equation  describing 
the  experimental  flow  properties  of  the  system  reveal  more  precise  outcomes  than 
previously  practiced  routines.  Resolved  viscosities  and  Poiseuille  flow  theory  velocity 
distributions  are  characterized  for  noble  gas  carriers  of  the  rubidium  vapor.  Applying 
Reynolds  Numbers  to  the  flow  experiments  provides  pressure  differential  boundaries 
that  are  employed  in  the  successful  rubidium  vapor  transfer  process.  Atomic  spec¬ 
troscopy  is  demonstrated  through  the  use  of  a  rubidium  Di  resonate  diode  laser  to 
record  an  absorption  spectrum  and  extract  alkali  vapor  number  densities  that  success¬ 
fully  propagated  through  a  12  cm  long  capillary  tube  with  a  500  fim  inner  diameter. 
Rubidium  number  densities  on  the  order  of  1.384  x  1012  cm”3,  4.615  x  1011  cm-3 
and  9.890  x  109  cm~3  were  recorded  for  interaction  path  lengths  of  0.05  cm,  0.15 
cm  and  7.0  cm,  respectively.  Number  densities  were  achieved  through  the  constant 
flow  and  diffusion  of  a  helium/rubidium  combination.  Sustaining  a  constant  pressure 
differential  between  80  and  150  Torr  across  the  capillary  tube  assisted  in  transferring 
the  rubidium  vapor  through  a  small  inner-diameter  hollow-core  tube. 
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DIFFUSION  OF  RUBIDIUM  VAPOR  THROUGH  HOLLOW-CORE  FIBERS 

FOR  GAS-PHASE  FIBER  LASERS 

I.  Introduction 

The  DPAL  and  Fiber  DPAL 

The  United  States  Air  Force  has  utilized  chemical  lasers,  specifically  the  chemical 
oxygen  iodine  laser  (COIL),  aboard  the  Airborne  Laser  (ABL)  and  the  Advanced 
Tactical  Laser  (ATL).  Solid  state  laser  systems  are  also  in  consideration  for  such 
projects,  but  they  introduce  thermal  management  issues.  Diode  pumped  alkali  lasers 
(DPAL)  have  recently  become  more  attractive  because  of  their  high  optical-to-optical 
efficiencies  and  low  thermal  loadings.  The  idea  behind  DPALs  is  to  turn  an  unphased 
diode  bars  into  a  powered  scaled  coherent  source.  These  alkali  vapor  systems  employ 
a  quasi-three  level  transition  scheme  where  a  narrow  banded  low  quality  beam  excites 
the  valence  electron  of  an  alkali  vapor  atom  from  the  2Si/2  (ground  state)  to  the  2P3/2 
excited  (D2)  state.  A  buffer  gas,  typically  ethane  or  methane,  is  used  to  collisionally 
relax  the  electron  to  the  2Pi/2  excited  (Dj)  state.  If  a  population  inversion  is  present 
on  the  Di  level,  the  system  will  lase  down  to  the  ground  state.  Due  to  the  alkali  vapor 
flowing  through  the  gain  cell,  thermo-optical  distortions  may  be  averted  providing  an 
output  beam  with  high  beam  quality. 

Increasing  the  flow  of  alkali  vapor  within  the  cell  has  been  exercised  through  the 
use  of  a  heat  pipe  to  introduce  more  atonr-photon  interactions.  A  heat  pipe  generates 
a  large  number  density  of  alkali  vapor  in  the  center  of  the  gain  cell.  As  the  temperature 
of  the  alkali  goes  up,  the  atoms  begin  to  flow  to  the  top  of  the  pipe.  Both  ends  of  the 
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Figure  1.  Heat  pipes  were  used  to  generate  large  number  densities  of 
alkali  in  DPAL  systems.  The  flow  of  the  alkali  allows  for  more  atomic- 
photon  interactions.  It  is  created  by  heating  the  center  of  the  pipe  and 
cooling  the  ends  with  flowing  chilled  water.  Buffer  gases  are  introduced 
to  collisionally  relax  the  excited  electrons  from  the  2P3/2  state  to  the 
2P 1/2  state. 


heat  pipe  are  sustained  around  room  temperature  conditions  or  slightly  lower  with 
the  aid  of  continuously  flowing  chilled  water,  which  attracts  the  hot  alkali  atoms. 
The  alkali  experiences  a  temperature  drop  as  it  approaches  the  ends  of  the  system 
and  will  eventually  begin  to  condense  and  relocate  to  the  bottom  of  the  cell  and  the 
process  is  repeated,  as  seen  in  Figure  1.  This  generates  a  flowing  system  within  the 
heat  pipe  allowing  for  more  alkali-photon  interactions  throughout  the  gain  medium. 
Constructing  open  and  closed  loop  continuously  flowing  alkali  systems  is  also  being 
investigated  to  generate  more  atomic-photon  interactions. 

Narrowing  the  diode  bars  is  needed  for  power  scaling.  Alkali  atomic  absorption 
profiles  are  extremely  narrow  requiring  the  diode  emission  bandwidth  to  also  be  nar¬ 
rowed  in  order  to  overlap  each  other  and  prevent  wasting  energy  and  promote  power 
scaling.  Creating  a  DPAL  system  within  a  fiber  would  create  an  abundant  number 
of  collisions  between  the  alkali  atoms  and  the  fiber  walls,  which  would  broaden  the 
atomic  spectral  width.  With  atomic-fiber  wall  collisional  broadening  in  consideration, 
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the  requirement  to  narrow  the  diode  emission  bandwidth  can  be  relaxed.  Utilizing 
the  waveguide  geometries  of  a  fiber  would  further  increase  the  absorption  of  the  diode 
because  of  the  large  length-to-diameter  ratios.  The  gain  medium  interaction  path  dis¬ 
tance  of  a  DPAL  system  could  be  increased  to  meters  rather  than  centimeters,  which 
is  currently  experienced  when  using  free  space  optics  due  to  the  multimode  nature 
and  required  intensities  of  the  diode.  Large  amounts  of  gain  over  longer  distances 
will  be  realized  from  the  large  surface  area  to  volume  ratios  characterized  by  fibers. 
Combining  these  characteristics  could  allow  a  shorter  operating  wavelength  region 
and  potentially  progress  to  the  development  of  a  true  hybrid  laser  system. 
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II.  Background 


Introduction 

This  work  examines  the  diffusion  of  rubidium  through  a  small  diameter  tube  alone 
and  in  the  presence  of  noble  gases.  A  diode  pumped  alkali  amplifier  laser  system  can 
be  created  if  enough  alkali  vapor  is  present  in  the  tube.  Previous  work  on  hollow-core 
fiber  lasers  will  be  mentioned  to  explain  the  different  types  of  lasers  that  have  been 
demonstrated  and  the  inner  diameters  that  have  been  used.  Understanding  the  fluid 
dynamics  taking  place  in  the  tube  will  avoid  chaotic  conditions  that  produce  shock 
waves  that  might  solidify  the  alkali  vapor.  Viscosity  measurements  are  necessary  to 
characterize  the  flow  within  the  tube  and  derivations  showing  how  to  experimentally 
extract  these  values  is  implemented.  Using  a  carrier  gas  to  assist  in  transferring  rubid¬ 
ium  vapor  through  the  capillary  tube  will  be  explored.  The  acquired  viscosity  values 
will  map  a  velocity  distribution  profile  determining  which  gas  has  the  highest  flow 
velocity  through  the  tube  while  avoiding  turbulent  flow  conditions.  An  understand¬ 
ing  of  atomic  rubidium  spectroscopy  is  developed  in  order  to  calculate  the  number 
density  of  alkali  present  in  the  experimental  system  after  the  transfer  process  has 
ocurred. 

Gas  Filled  Hollow-Core  Fiber  Lasers 

Techniques  for  generating  high  number  densities  of  atomic  or  molecular  gases  in 
a  fiber  are  recently  in  high  demand.  The  practice  of  filling  hollow-core  fibers  with 
gases  for  lasing  purposes  is  unique  because  it  conjoins  the  techniques  of  light-atom 
interactions  with  the  waveguide  properties  of  a  fiber.  The  waveguide  geometry  of  the 
fiber  constrains  the  light  to  propagate  through  the  gain  medium,  such  as  alkali  vapor, 
over  long  interaction  lengths  with  little  loss  due  to  scatter  and  transmission  through 
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the  fiber.  Gas  filled  fiber  lasers  also  contribute  to  the  concept  associated  with  DPALs 
of  generating  a  coherent  beam  from  an  incoherent  source.  The  cardinal  applications 
for  gas  filled  fibers  include  atomic  spectroscopy,  quantum  interference,  and  slow  light 
experiments. 

Slepkov  et  al.  [7]  generated  high  rubidium  number  densities  through  diffusion 
inside  a  4.5  cm  long,  bare  hollow-core  photonic  band-gap  fiber  with  a  core-diameter 
of  6  nrm.  Number  densities  on  the  order  of  108  cm~ 3  were  achieved  in  the  fiber 
using  light-induced  atomic  desorption  (LIAD).  L1AD  is  a  technique  which  uses  a 
weak  resonate  laser  to  remove  measurable  atomic  densities  from  an  alkali-exposed 
glass  surface.  Diffusion  processes  will  be  discussed  in  this  thesis  in  attempt  to  gain 
a  number  density  on  the  order  required  to  operate  a  DPAL  system,  which  is  on  the 
order  of  1013  cm~3.  Smith  [8]  was  the  first  to  report  a  helium-neon  laser  system 
within  a  small  inner- diameter  (ID)  tube  in  1971.  Using  a  20  cm  long  glass  capillary 
tube  with  a  430  / im  ID,  Smith  used  a  633  nm  laser  source  to  create  a  helium-neon 
fiber  laser.  Jin  et  al.  [3]  exhibited  discharges  of  argon,  helium  and  carbon  dioxide  in 
hollow-core  fibers  with  IDs  of  250  /im,  150  gm,  and  50  fmi  and  a  length  up  to  13.7 
cm.  Jones  et  al.  [4]  demonstrated  an  optically  pumped  gas  laser  in  a  hollow-core 
photonic  crystal  fiber  using  acetylene  (C2H2).  A  pulsed  source  providing  5  —  6  ns 
pulses  of  approximately  5  mJ  of  energy  per  pulse  was  directed  into  a  1.65  m  long 
single-cell  optical  fiber  containing  C2H2.  A  population  inversion  was  confirmed  as  two 
rotational  branch  transitions  ((i?7)  and  (-P9))  were  observed  near  3.12  and  3.16  fim. 
A  slope  efficiency  of  a  few  percent  was  achieved  even  though  the  fiber  coupling  loses 
were  around  98%.  Hostutler  et  al.  [2]  demonstrated  a  rubidium  vapor  laser  with  a 
master  oscillator  power  amplifier.  An  amplification  of  5.9  d.B  and  7.9  dB  was  achieved 
using  a  1.5  cm  and  2.0  cm  cell,  respectively.  The  master  oscillator  performed  at  a 
low  intensity  to  avoid  saturating  the  sample.  A  14%  maximum  extraction  efficiency 
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was  determined  fr  the  alkali  amplifier. 


Capillary  Tube  Fluid  Dynamics 

As  various  gases  propagate  through  the  capillary  tube  (a  surrogate  for  a  hollow- 
core  fiber)  a  series  of  fluid  dynamical  phenomena  occur.  Pressure  measurements 
reveal  flow  characteristics  such  as  diffusion  rates,  gas  viscosity,  flow  velocity,  and  flow 
velocity  profiles.  Understanding  the  fluid  dynamics  will  assist  in  the  selection  of  which 
gas  to  choose  to  efficiently  aid  in  transferring  a  number  density  of  alkali  vapor  through 
a  small  inner  diameter  capillary  tube.  Experimentally  determining  the  gas  viscosity 
requires  the  ability  to  model  the  gas  diffusion  rate  as  it  flows  through  the  tube.  The 
following  derivations  discuss  how  the  viscosities  can  be  estimated  using  variations  of 
the  ideal  gas  law,  the  geometry  of  the  capillary  tube,  and  the  pressure  of  the  chambers 
on  both  ends  of  the  tube.  The  first  avenue  explores  a  common  technique  that  assumes 
a  large  pressure  on  the  gas  inlet  side  of  the  tube  and  a  negligible  pressure  in  the  aft 
chamber.  The  second  derivation  will  include  the  previously  neglected  back  chamber 
pressure  and  outline  all  of  the  physical  properties  occurring  in  the  flow  analysis.  Once 
the  viscosities  are  ascertained,  the  Reynolds  Number  will  be  introduced  to  determine 
what  type  of  flow  is  being  demonstrated  in  the  tube. 

The  two  most  common  types  of  flow  for  a  cylindrical  tube,  like  the  one  used 
in  this  experiment,  are  turbulent  flow  and  laminar  flow.  When  turbulent  flow  is 
present,  the  gas  travels  down  the  pipe  in  a  chaotic  and  an  unpredictable  manner 
with  an  uneven  velocity  profile  distribution.  It  is  difficult  to  observe  and  understand 
the  dynamics  of  a  turbulent  flow  scenario.  Laminar  flow  theory,  also  referred  to  as 
Poiseuille  flow  theory,  is  predictable  and  provides  an  easier  analysis  of  the  velocity 
profile  distribution.  When  laminar  flow  is  present,  the  gas  near  the  center  of  the 
capillary  tube  travels  with  a  higher  velocity  than  on  the  outer  edges  of  the  tube, 
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represented  by  v(r)  in  Figure  2.  Because  the  wall  of  the  tube  is  stationary,  it  applies 
a  small  frictional  force  on  the  flowing  gas.  The  velocity  of  the  gas  at  the  wall  of  the 
tube  is  zero  during  laminar  flow.  The  frictional  force  experienced  on  the  adjacent 
laminar  sheets  as  they  flow  by  can  be  described  by: 

*  =  =  (1> 

where  rj  is  the  viscosity  of  the  gas,  A  is  the  area  of  contact  between  the  laminar  sheets, 
l  is  the  length  of  the  tube,  v  is  the  flow  velocity,  and  r  is  the  distance  from  the  center 
of  the  capillary  tube,  making  ^  the  velocity  gradient  through  the  cross-section  of  the 
tube.  The  force  driving  the  fluid  through  the  tube  is: 

fd=(Pi~P2)irr2,  (2) 

where  p\  is  the  fore  pressure  and  P2  is  the  back  pressure.  When  steady  flow  is  observed 
through  the  tube,  the  forces  can  be  described  as  fully  developed  and  counteract  one 
another  when  summed  together: 


fr  +  fd  —  0 


(3) 


dv 

r)(2nrl)—  +  (pi  -  p2)7rr2  =  0. 
dr 


(4) 


The  fact  that  the  opposing  forces  sum  to  zero  is  a  consequence  of  the  law  of  momentum 
conservation  of  the  gas  propagating  through  the  capillary  tube.  By  solving  for  the 
change  in  the  velocity  of  the  gas  and  reducing  the  equation  the  result  is: 


dv  =  "fo-P’W. 
2/7 1 


(5) 
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(a) 


< 


(b) 


a  r &a  =  Zsrdr 


Figure  2.  (a)  A  pipe  experiencing  a  pressure  gradient  ( pi  —  p2)  and 

under  laminar  flow  (Poiseuille  flow)  has  a  parabolic  velocity  profile.  The 
velocity  in  the  center  of  the  tube  is  the  highest  and  the  velocity  near  the 
walls  is  the  lowest.  As  the  parabolic  profile  travels  through  the  capillary 
tube  the  gas  velocity  remains  constant  along  the  length  of  the  tube,  (b) 
Cross-section  of  the  capillary  tube  showing  geometrical  dimensions. 


In  order  to  know  the  fluid  velocity  as  a  function  of  distance  from  the  center  of  the 
capillary  tube  (he.  the  velocity  profile  distribution)  the  above  equation  is  integrated 
from  the  edge  of  the  pipe  to  the  center  of  the  pipe  (v  =  0  at  R  —  r,  where  R  is  the 
radius  of  the  capillary  tube): 


dl,  =  -(fi-p2)  rrdr 


2rjl 


(6) 


lR 


v(r )  = 


(Pi  -  P2)  ,  D2  2n 


477/ 


-(R2  -r2). 


(7) 


equation  (7)  describes  the  gas  velocity  profile  propagating  down  the  center  of  the  tube 
as  a  paraboloid  when  laminar  flow  is  present,  shown  in  Figure  2a.  As  the  parabolic 
profile  travels  through  the  capillary  tube  the  gas  velocity  remains  constant  along  the 
length  of  the  tube. 

The  flow  rate  through  a  tube  connecting  two  chambers  under  a  pressure  gradi¬ 
ent  can  be  found  by  understanding  the  geometry  of  the  capillary.  By  knowing  the 
dimensions  of  the  tube,  the  volumetric  flow  rate  can  be  defined  as: 


dV 

dt 


rR 


2nvrdr, 


(8) 


where  V  is  the  volume  containing  the  fore  pressure.  A  cross-section  of  a  cylindrical 
tube  is  shown  in  Figure  2b.  By  substituting  equation  (7)  into  equation  (8),  the  change 
of  volume  of  gas  per  time  becomes: 


dV 

dt 


(27r)A  p 
A'ql 


r2)dr, 


(9) 


where  A p  is  the  pressure  gradient  across  the  tube,  (pi  — p2)-  The  volumetric  flow  rate 
can  be  reduced  to: 


dV  _  vr(A p)R4 
dt  8r/l 
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If  the  average  pressure  of  the  system  is  defined  as: 


P 


1 

2 


(P1+P2), 


(11) 


and  if  the  gas  is  assumed  to  be  ideal,  the  number  of  moles  of  gas  flowing  through  the 
tube  per  unit  time  is: 


dn  P  dV 
dt  ksT  dt  ’ 


(12) 


where  n  is  the  number  of  moles  of  flowing  gas,  ks  is  the  Boltzmann  constant,  and  T  is 
the  absolute  temperature  of  the  system.  By  substituting  in  the  volumetric  flow  rate 
and  multiplying  by  one  using  the  average  pressure  with  the  motive  of  cancellation, 
the  result  becomes: 


dn  _  P  nApR 4  (pi  +  p2) 
dt  ksT  8r)l  2  P 

Once  the  gas  is  added  to  the  fore  pressure  chamber,  it  can  be  assumed  that  the  fore 
pressure  is  much  greater  than  the  back  pressure  {p\  3>  P2)  and  the  flow  of  gas  per 
unit  time  is  found  through: 


dn  tv  p\RA 

dt  16  ksTijl 


(14) 


The  ideal  gas  law  describes  the  change  of  the  fore  pressure  per  unit  time  by: 


dpi  —ksT  dn 
dt  V  dt' 


and  by  inserting  equation  (14)  into  equation  (15),  the  statement  becomes: 


dpi  —  7ip\RA 
dt  16  J]IV 


The  diffusion  rate  through  the  capillary  tube  can  be  described  by  the  separation  of 
variables  in  equation  (16)  and  by  integrating  the  fore  pressure  from  the  initial  pressure 
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to  the  instantaneous  fore  pressure  and  by  integrating  over  all  experimental  time.  The 
expression  can  be  written  as: 


fPldp  i  itR4  f* 

JPi  pi  iQvivJo  ’ 

which  can  be  simplihed  to  the  reciprocal  of  the  fore  pressure  as: 

1  _  ttR4  1 

Pi(t )  16  7]lVt  + p.i' 


(17) 


(18) 


where  Pi  is  the  initial  pressure  of  the  system  and  t  is  time.  If  the  geometrical  dimen¬ 
sions  of  the  capillary  tube  are  known,  the  time  dependant  derivative  of  equation  (18) 
provides  an  expression  that  aids  in  finding  the  viscosity  of  the  gas  flowing  through 
the  tube  experimentally: 


d(±)  _  t tR4 
dt  16plV ' 


(19) 


Figure  3  shows  the  reciprocal  of  simulated  experimental  fore  pressure  data  plotted 
versus  time.  The  viscosity  of  the  gas  can  be  established  by  fitting  the  slope  of  the 
initial  portion  of  the  inverse  fore  pressure  data.  The  previous  derivations  demonstrate 
a  technique  for  determining  the  viscosity  which  will  assist  in  characterizing  the  type 
of  flow  in  the  capillary  tube  [5,6]. 

Rather  than  exclusively  fitting  the  slope  of  the  inverse  pressure,  more  accurate 
viscosity  results  can  be  extracted  by  solving  the  ordinary  differential  equation  that 
describes  the  flow  through  the  capillary  tube  from  one  chamber  to  the  other.  By 
using  the  ideal  gas  law  in  equation  (15)  and  solving  for  ^  and  inserting  the  result 
into  equation  (13),  the  expression  becomes: 


dpi  _  n(pi  -p2)R4  (pi  +P2 ) 
dt  8  rjVl  2 


(20) 
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Pl~  P2 


1/Pl 


Time 


Figure  3.  Simulated  data  of  the  inverse  fore  pressure  plotted  versus  time. 
The  initial  slope  of  the  curve,  before  the  roll-off  where  the  fore  pressure 
is  approximately  equivalent  to  the  back  pressure,  provides  an  analytical 
expression  that  can  be  solved  for  the  flowing  gas  viscosity. 


The  assumption  where  p\  P2  will  be  ignored  in  order  to  consider  all  the  physical 
properties  of  the  system.  The  pressure  that  describes  the  total  pressure  of  the  system 
can  be  defined  by: 


(21) 


Pi)  =  Pi  +  P2- 


Substituting  the  back  pressure,  p2,  into  equation  (20),  the  ordinary  differential  equa¬ 
tion  (ODE)  describing  the  fore  pressure  as  a  function  of  time  becomes: 


Pi(t) 


tt(Pi^)  -  (Po  ~  Pi(t)))R4 Pi(t)  +  (p0-Pi(t)) 


(22) 


8  rjVl 


2 


The  solution  to  the  ODE  is: 


(23) 
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The  objective  is  to  fit  the  entire  inverse  fore  pressure  curve  which  is: 


1 

Pi(t) 


2  exp 


P0irR4t 

8r]Vl 


Pofi  +  ^pf'lfr 


(24) 


Using  the  MATLAB  curve  fitting  function  (which  utilizes  the  Least-Squares,  Levenberg- 
Marquardt,  and  Gauss-Newton  methods)  the  fit  parameters  are  the  total  pressure  of 
the  system,  p0,  and  T,  where: 


r 


7 lR4 

8pVl 


(25) 


This  makes  the  fit  function  (with  only  two  fit  parameters): 


1 

Pi(t) 


Po 


2exp(poTt) 

^1  +  exp(p0Tt) 


(26) 


where  po  is  expressed  in  Torr  and  T  is  expressed  in  ( Torr  —  Sec )_1. 

The  previous  equations  demonstrate  the  derivations  leading  to  discovering  the 
viscosity  of  the  flowing  gas  as  a  paraboloid  through  the  center  of  the  capillary  tube 
under  the  conditions  of  laminar  flow  theory.  In  order  to  determine  whether  the  flow 
is  turbulent  or  laminar,  the  Reynolds  Number  must  be  calculated  and  analyzed  for 
the  specific  gas  under  investigation.  The  Reynolds  Number  can  be  found  by: 


Re  = 


(27) 


where  p  is  the  density  of  the  gas,  v  is  the  average  velocity  of  the  gas  through  the 
tube,  d  is  the  diameter  of  the  tube,  and  p  is  the  viscosity  of  the  gas  determined 
experimentally  from  equation  (19)  or  equation  (26).  The  average  velocity  of  the  gas 
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is  defined  as: 


v  = 


(Pi  ~P2)R2 


8r]l 


(28) 


If  the  Reynolds  Number  is  less  than  approximately  30,  the  Poiscuille  flow  theory 
always  provides  an  accurate  description  of  the  dynamic  flow  [10,11].  At  higher  values 
of  the  Reynolds  Number,  the  Poiseuillc  flow  theory  applies  only  after  some  distance 
down  the  tube  as  it  is  unlikely  the  gas  enters  the  tube  with  the  appropriate  parabolic 
velocity  profile  defined  in  equation  (7).  If  the  Reynolds  Number  is  outside  of  the 
Poiseuillc  flow  regime,  the  distance  down  the  tube  before  it  follows  the  parabolic 
velocity  profile  is  calculated  by: 


X 


(. Re)d 
30 


(29) 


If  the  pressure  gradient  is  equivalent  for  each  gas  under  inspection  and  the  radius 
and  length  of  the  tube  are  not  changed,  the  maximum  velocity  of  the  gases  is  solely 
dependent  on  the  viscosity  of  each  gas.  Therefore,  it  is  critical  to  distinguish  which 
type  of  flow  is  occurring  in  the  system  so  that  the  experimental  viscosity  results  can 
be  confirmed  or  discarded.  The  maximum  velocity  of  the  gas  flowing  through  the 
tube  is: 

(pi  ~P-2)R2 

vmax  = - - 

4r/ 

and  can  assist  in  the  selection  of  which  gas  to  choose  to  efficiently  aid  in  transferring 
a  number  density  of  alkali  vapor  through  a  small  inner-diameter  capillary  tube. 


Multi-isotopic  Rubidium  Spectroscopy 

It  is  imperative  to  recognize  the  fundamental  quantum  principals  when  investigat¬ 
ing  an  atomic  rubidium  spectrum.  The  core  objective  pertaining  to  this  experiment 
is  to  determine  the  rubidium  vapor  number  density  in  the  capillary  tube  confirming 
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a  sufficient  amount  of  alkali  is  present  within  the  tube  to  generate  a  diode  pumped 
alkali  amplifier.  Rubidium  has  two  naturally  occurring  isotopes,  85 Rb  and  8'Rb  with 
abundances  72.17%  and  27.83%,  respectively.  The  fine  structure  of  the  rubidium 
atom  is  defined  by  coupling  the  orbital  angular  momentum,  L ,  with  the  spin  angular 
momentum,  S,  to  provide  the  total  angular  momentum  of  each  fine  structure  energy 
level,  J  =  L  +  S.  The  energy  difference  between  the  fine  structure  excited  states  for 
both  isotopes  is: 

Efs  =  E(2P3/2 )  -  E(2P1/2)  =  7.123  020  683  THz  ±  60  kHz.  (31) 

When  the  total  angular  momentum  is  coupled  with  the  nuclear  spin,  /  (/  =  5/2 
for  85Rb  and  /  =  3/2  for  8'Rb),  the  hyperfme  structure  of  the  rubidium  doublet  is 
resolved.  The  hyperfine  structure  energy  levels  follow  the  quantum  selection  rules  as: 

\J-I\  <  F  <  J  +  I,  (32) 


where  F  is  the  total  atomic  angular  momentum.  The  hyperfine  structure  energy 
splitting  is  defined  by: 


Ehfs(F)  —  Efs  +  ^  +  K 


(3/2)x(x  +  1)  ~  21(1  +  1)J(J  +  1) 
4/(2/-  1)J(2J-  1) 


(33) 


where 


and 


0) 

IJ 


magnetic  dipole  constant, 


k  =  eQ4>jj( 0)  =  electric  quadrupole  constant, 


(34) 

(35) 


y  =  F(F  +  1)  —  /(/  +  1)  —  J(J  +  1)  =  magnetic  octupole  constant.  (36) 
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It  must  be  noted  that  the  terms  in  equation  (33)  containing  k  exclusively  apply  to 
the  energy  states  with  a  total  angular  momentum  of  J  =  3/2  and  do  not  apply  to  the 
Di  transition.  The  hyperhne  structure  energy  level  splitting  for  both  of  the  rubidium 
isotopes  are  shown  for  the  Di  line  in  Figure  4  and  Figure  5.  As  a  Di  laser  source 
scans  across  the  resonate  transition  peaks  of  the  multi-isotopic  rubidium  vapor,  near 
infrared  photo-detectors  can  be  used  to  capture  the  atomic  absorption  spectrum.  If 
the  probe  intensity  is  less  than  the  saturation  intensity  ( I/Isat  <  1)  and  the  sample 
is  not  being  bleached,  the  absorbance  across  the  test  cell  can  be  described  through 
Beer’s  Law  as: 

Aa  =  -  ln(/t//0)  =  a(u,  T)lc  =  a(u,  T)N(T)lc,  (37) 


where  It  is  the  transmitted  laser  intensity  through  the  test  cell,  J0  is  the  intensity  of 
the  laser  entering  the  test  cell  (background  intensity),  a  is  the  absorption  coefficient, 
v  is  the  laser  frequency,  T  is  the  temperature  of  the  cell,  lc  is  the  length  of  the  cell,  a  is 
the  cross-section  of  absorption  and  N  is  the  rubidium  number  density  in  the  sample. 
The  //Jo  data  will  be  extracted  from  the  experiment  to  determine  the  absorbance.  If 
the  length  of  the  cell  is  known  and  the  absorption  cross-section  can  be  determined, 
the  number  density  of  the  system  can  be  ascertained.  The  cross-section  of  absorption 
is  found  through: 


2vr2 

<7F"F'(v,T)  =  - — —ii2F„F,h/g(u  —  uF"Fi ) , 

oCoflC 


(38) 


where  eo  is  the  vacuum  permittivity,  h  is  Planck’s  constant,  c  is  the  speed  of  light  in 
a  vacuum,  fiF”F>  is  the  dipole  moment  and  g(u  —  uF”Fr)  is  the  absorption  lineshape 
at  a  resonate  frequency  transition  uF" F>  or  //0-  The  dipole  moment  can  be  found  by 
applying  the  dipole  moment  operator  on  the  hyperhne  transition  under  investigation: 


—  SF"F'\{J"\\er\\J')\2  —  SF"Fie2D2j/ij,. 


(39) 
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35Rb 


F=  3 


361.58  MHz 


F=  2 


F=  3 


3.035  732  439  GHz 


F=  2 


Figure  4.  The  Grotrian  diagram  for  the  85Rb  Di  line  is  shown  in  including 
the  hyperflne  structure  and  energy  level  splitting  frequencies  [9]. 


17 


a7Rb 
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52S 


F=  2 


816.656  MHz 


F=  1 


F=  2 


6.834  682  GHz 


F=  1 


Figure  5.  The  Grotrian  diagram  for  the  8'Rb  Dx  line  is  shown  in  including 
the  hyperflne  structure  and  energy  level  splitting  frequencies  [9]. 
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The  relative  line  strength,  Sp"F is  resolved  using  Wigner’s  6-j  matrix  following: 


S 


F"F' 


(2F'  +  1)(2J"  +  1) 


J" 

F' 


J' 

F" 


(40) 


Relative  line  strengths  with  their  associated  upper  (. F ')  and  lower  (F")  hyperfine 
states  are  defined  for  each  isotope  in  Table  1.  The  dipole  moment  is  used  to  ascer- 


Table  1.  Hyperfine  transitions  following  quantum  selection  rules  for  the 
total  atomic  angular  momentum  (equation  (32))  for  each  isotope  are 
associated  with  a  line  strength  ( Sf"f ')  and  a  detuning  frequency  (v  —  vo) 

[9] 


S f" f'  (unitless) 

v-v0 

85Rb 

F"  =  2,  F' 

=  3 

0.778 

1.921  502  GHz  ±  71  kHz 

F"  =  2,  F' 

=  2 

0.222 

1.559  920  GHz  ±  99  kHz 

F"  =  3,  F' 

=  3 

0.444 

-1.114  229  GHz  ±  71  kHz 

F"  =  3,  F' 

=  2 

0.556 

-1.475  811  GHz  ±99  kHz 

87Rb 

F"  =  1,  F' 

=  2 

0.833 

4.577  922  GHz  ±  11  kHz 

F"  =  1,  F' 

=  1 

0.167 

3.761  266  GHz  ±  19  kHz 

F"  =  2,  F' 

=  2 

0.500 

-2.256  759  GHz  ±  11  kHz 

F"  =  2,  F' 

=  1 

0.500 

-3.073  415  GHz  ±  19  kHz 

tain  the  radiative  lifetime  of  the  each  energy  level  of  the  rubidium  atom,  but,  the 
degeneracy  must  be  considered  to  differentiate  the  levels.  The  arithmetic  to  calcu¬ 
late  the  dipole  moment  can  be  quite  involved,  but  can  become  relatively  trivial  when 
the  dipole  moment  is  substituted  for  the  radiative  lifetime.  In  terms  of  the  dipole 
moment,  the  radiative  lifetime  is: 


1 

Tr 


gj" 

S7re0hc 3  gj> 


e2D 


2 

J"J'i 


(41) 


where  ojq  is  the  resonate  angular  fequency  (cuo  =  27t u0),  h  is  Planck’s  constant  divided 
by  27 r,  and  gj  is  the  degeneracy  of  the  energy  level,  J .  This  simplifies  the  dipole 
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moment  to  become: 


2  3eo he3  /  2J7  +  1  \  1 

and  the  cross-section  of  absorption  reduces  to: 


(42) 


(  ^r\  c2  v  f  2  J’  +  1  ^  1 

<jF"f'{v,T)  -  —  -3  0  T//  ,  )SF"F’—g(v  -  v 0). 


87t  vS  \2  J"  +  1 


Tr 


(43) 


The  Boltzmann  distribution  considers  the  probability  for  a  single  rubidium  atom 
to  be  excited  to  one  of  the  hyperhne  structure  energy  states  that  are  shown  in  Figure 
4  and  Figure  5.  These  probabilities  are  based  on  the  atomic  structure  and  degenera¬ 
cies  of  the  atom  as  well  as  the  lower  hyperhne  energy  level  of  the  transition.  The 
Boltzmann  distribution  can  be  expressed  by: 


,/f" 


(2  F"  +  l)e~EF"/kBT 

£<2F"  +  X)e~Ep"  1 

F" 


(44) 


where  EFn  is  the  energy  difference  of  the  52S'i/2  hyperhne  quantum  splitting  and 
F"  is  the  total  atomic  angular  momentum  of  the  lower  transition  state.  The  rela¬ 
tive  amplitude  for  absorbance,  —  ln(/t//o),  for  each  resonate  transition  is  simply  the 
product  of  the  line  strength,  SF» F> ,  the  Boltzmann  distribution,  fFn:  and  the  natural 
isotopic  abundance.  The  relative  amplitudes  of  each  transition  are  shown  in  Figure 
6  at  room  temperature  and  will  assist  in  fitting  the  acquired  rubidium  spectrum.  In¬ 
serting  the  distribution  across  all  of  the  possible  resonate  energy  states  reduces  the 
total  cross-section  of  absorption  to: 


a(V,T)  =  £  fF„S 

F"  ,F' 


1  c2  V 


(45) 


The  lineshape  of  the  atomic  transitions  can  be  expressed  by  the  Voigt  prohle 
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Figure  6.  Relative  frequency  locations  (v—Vq)  and  relative  amplitudes  for 
each  hyperfine  transition  at  room  temperature.  The  relative  amplitude 
is  the  product  of  the  line  strength,  Sf"f'i  the  Boltzmann  distribution, 
fp"t  and  the  natural  isotopic  abundance. 


approximation: 


9v{v 


Vo)qh(v  ~  v'0)dv'0, 


(46) 


which  contains  a  homogeneous  and  an  inhomogeneous  part  [1].  The  homogeneous 
part  contains  information  about  the  natural  linewidth  of  the  atom  and  any  pressure 
broadening  experienced  in  the  system: 


9 ~  Vo) 


(A^h/2)2  +  (v  -  vaf' 


(47) 


where  the  full-width  at  half  maximum  (FWHM)  is  illustrated  by  a  Lorentzian  line- 
shape  denoted  by: 


T 


,7* 


T) 


1/2 


T, 


Pi. 


data 


(48) 


27T  \Tr/ 

The  pressure  broadening  factor  contains  the  broadening  coefficient,  y*,  the  tempera¬ 
ture  at  which  the  broadening  coefficients  were  calculated  in  order  to  ration  the  amount 


21 


of  broadening,  Tdata,  and  the  pressure  of  the  system,  Pt. 

The  inhomogeneous  portion  of  the  Voigt  profile  approximation  can  be  determined 
from: 


9i{v“Vo)  =  -£ 2~i^eXV 


-  4  In  2 


v-vQ 


Az/, 


G 


(49) 


which  is  represented  by  a  Gaussian  lineshape  and  the  FWHM  is  presented  as: 


A  vG  =  v0 


8k bT  In  2\ 
Me 2  )’ 


(50) 


where  M  is  the  atomic  mass  of  rubidium  (85Rb=  1.409-25  kg,  8'Rb=  1.4432”25  kg) 
[9].  The  Voigt  profile  is  the  convolution  of  the  Lorentzian  and  Gaussian  lineshapes. 
With  the  system  lineshape,  the  absorption  cross-section  in  equation  (45)  can  be  solved 
using  known  experimental  and  spectroscopic  information  [1],  The  low  pressure  cross- 
section  at  line  center,  <r0  =  1.081  257  448  975  x  10-9  ±  54  x  10-21  cm2,  will  be  used 
to  assist  in  resolving  a  number  density.  Experimental  procedures  will  provide  the 
amount  of  absorbance,  A\,  and  the  alkali  vapor  number  density  generated  in  each 
chamber  can  be  extracted  using  equation  (37).  If  a  number  density  can  be  extracted 
from  the  alkali  chambers  at  each  end  of  the  capillary  tube,  a  steady  number  density 
gradient  can  be  assumed  across  the  tube  to  determine  whether  an  alkali  amplifier  is 
feasible. 
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III.  Experimental  Setup 


Flow  Rate  Experiments 

Understanding  the  fluid  dynamics  of  various  gases  as  they  propagate  through  the 
capillary  tube  will  assist  in  choosing  the  most  efficient  carrier  for  the  rubidium  vapor. 
By  determining  the  Reynolds  Numbers  and  confirming  the  flow  is  characterized  by 
Poiseuille  flow  theory,  shock  waves  and  chaotic  flow  behavior  common  during  turbu¬ 
lent  flow  conditions  are  averted.  Shock  waves  have  the  potential  of  solidifying  the 
rubidium  if  large  enough  pressures  are  present  within  the  tube.  Before  the  rubidium 
was  added  to  the  experiment,  flow  rates  were  observed  from  four  of  the  noble  gases. 

Two  stainless  steel  cubic  chambers  with  an  internal  volume  of  158.79  cm3  were 
joined  together  by  a  12  cm  long,  quarter  inch  diameter  Polymicro  glass  rod  with  a 
500  firn  hole  through  the  center  (#TQNC5006350/200),  displayed  in  Figure  7.  The 
rod  was  conjoined  to  the  steel  cubes  using  an  Ultra- Torr-to-conflat  adapter  so  both 
ends  of  the  rod  were  located  in  the  center  of  the  chambers.  A  lecture  bottle  gas  rack 
was  fabricated  in  a  fashion  so  helium,  argon,  krypton,  and  xenon  could  be  added 
at  separate  times  by  turning  the  appropriate  shut-off  valves.  The  main  bus  of  the 
noble  gases  was  attached  to  the  top  of  one  of  the  cubes  acting  as  the  fore  pressure 
chamber  and  the  gases  could  be  added  gradually  using  a  needle  valve.  The  top  of 
the  back  pressure  cube  was  affixed  to  a  Varian  turbo  pump  that  achieved  pressures 
as  low  as  60  nTorr.  The  vent  of  the  turbo  pump  was  connected  to  a  roughing  pump 
that  disposed  of  the  gases  pumped  out  of  the  system.  Two  1000  Torr  MKS  Baratron 
(type  631)  pressure  sensors  were  attached  to  the  top  of  each  cube  and  were  wired 
to  a  PC  using  a  National  Instruments  (. NI  —  USB  —  6221)  DAQ,  which  recorded 
the  pressure  of  each  cube  simultaneously  using  Lab  View.  Using  shut-off  valves,  each 
chamber  could  be  isolated  from  the  gas  inflow  and  turbo  pump,  respectively. 
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Figure  7.  Fluid  dynamics  experimental  setup  with  a  12  cm  long,  quarter 
inch  Polymicro  glass  rod  with  a  500  fim  hole  through  the  center.  Each 
end  of  the  tube  is  located  in  the  center  of  both  chambers.  Two  1000  Torr 
MKS  Baratron  (type  631)  pressure  sensors  were  attached  to  the  top  of 
each  chamber.  A  gas  line  capable  of  dispensing  four  noble  gases  was 
attached  to  the  top  of  the  fore  pressure  cube  and  a  Varian  turbo  pump 
was  attached  to  the  top  of  the  back  pressure  cube. 


A  leak  check  was  conducted  by  connecting  the  system  to  a  helium  leak  detector 
and  dispensing  helium  on  all  of  the  connections  of  the  system  until  there  were  no 
apparent  leaks.  The  back  pressure  chamber  was  pumped  down  to  low  pressures  then 
was  shut  off  from  the  turbo  pump.  An  initial  pressure  of  each  gas  was  added  to  the 
fore  pressure  chamber  at  different  times.  The  gas  diffused  through  the  glass  rod  filling 
the  back  pressure  cube  until  both  were  at  equilibrium.  Pressures  from  both  volumes 
were  recorded  during  the  diffusion  process. 

Transferring  Rubidium  Vapor  Through  a  Capillary  Tube 

Nobel  gases  were  introduced  into  to  the  fore  pressure  cube  using  the  gas  rack  de¬ 
scribed  previously.  The  gases  were  proposed  as  potential  carriers  to  aid  in  transferring 
rubidium  vapor  atoms  through  the  tube  and  into  the  pump  side  chamber.  Generat¬ 
ing  a  large  alkali  vapor  number  density  for  the  diffusion  process  required  heating  the 
system  to  temperatures  ranging  from  60  C  to  150  C.  Heating  the  cubes  and  capillary 
tube  as  uniformly  as  possible  was  essential  to  ensure  no  rubidium  accumulated  in  cold 
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areas  of  the  system.  A  steel  housing  was  placed  around  the  system  with  six  radiating 
heaters  affixed  in  a  formation  to  adequately  increase  the  temperature  of  the  chambers 
and  the  capillary  tube  in  a  uniform  manner.  Both  of  the  chambers  were  placed  on 
half-inch  thick  metal  plates.  Each  plate  contained  four  equally  spaced  holes  burrowed 
horizontally  with  diameters  just  large  enough  for  current  driven  heater  elements  to  be 
inserted.  The  heater  elements  and  resistance  temperature  detectors  (RTD)  placed  on 
the  bottom  of  each  plate  were  connected  to  a  modified  Watlow  heater  controller  that 
prevented  increasing  and  decreasing  temperature  ramping.  Thermocouples  probed 
the  interior  of  the  housing  allowing  an  accurate  temperature  reading  of  the  radiation 
being  projected  onto  the  exterior  of  the  cubes  and  capillary  tube.  Calcium  fluoride 
windows  were  positioned  in  six  locations  where  diode  laser  beams  entered  or  exited 
the  system.  A  cold  trap  filled  with  liquid  nitrogen  was  placed  at  the  exit  of  the  back 
pressure  chamber  as  a  safety  net  to  solidify  and  condense  any  rubidium  atoms  before 
reaching  the  turbo  pump.  The  vent  of  the  turbo  pump  was  connected  to  a  roughing 
pump  to  dispose  of  the  substances  extracted  from  the  experimental  apparatus. 

Once  the  system  was  heated  to  create  rubidium  vapor,  three  techniques  were 
demonstrated  with  the  main  goal  of  transferring  the  alkali  through  the  capillary  tube 
and  into  the  pump  side  chamber.  The  Erst  procedure  was  to  continuously  pump  on 
the  back  pressure  chamber  with  the  turbo  pump  while  the  fore  pressure  chamber 
exclusively  contained  rubidium  vapor.  The  second  routine  exploited  basic  diffusion 
principals  through  the  inclusion  of  some  noble  gases  (He,  Ar,  Kr  and  Xe).  After 
pumping  the  system  down  to  low  pressures  (■ mTorr ),  the  gases  were  independently 
added  to  the  fore  pressure  chamber  combining  with  the  alkali  and  diffusing  to  the 
pump  side  cube.  The  third  approach  consolidated  the  two  previously  discussed  con¬ 
cepts.  A  noble  gas  was  continually  introduced  to  the  fore  pressure  chamber  containing 
the  alkali  vapor  while  the  valve  to  the  turbo  pump  was  slightly  opened  on  the  pump 
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side.  Keeping  both  input  gas  and  output  pump  gateways  barely  open  maintained  a 
deliberate  pressure  gradient  across  the  capillary  tube  that  was  sustained  by  adjusting 
the  needle  valves  as  necessary.  A  mixture  of  the  alkali  and  the  noble  gas  flowed  con¬ 
tinuously  through  the  tube  and  into  the  pump  side  chamber.  Each  of  these  avenues 
were  explored  and  each  achieved  their  own  rates  of  success. 

Detecting  a  Rubidium  Spectrum 

A  diode  laser  was  employed  to  scan  across  the  rubidium  Di  resonate  peaks  (ap¬ 
proximately  795  mm)  and  was  driven  by  a  Newport  diode  driver  (560B)  and  a  Thor- 
Labs  temperature  controller  (TED200).  The  diode  output  80  mW  of  power  with  a 
spectral  linewidth  FWHM  of  2  MHz  and  has  a  spot  size  of  approximately  0.5  cm. 
The  saturation  intensity  for  rubidium  is  4.4876  ±  0.0031 mW/cm2,  which  will  intro¬ 
duce  complications  when  fitting  the  spectrum  and  will  be  discussed  later.  A  signal 
generator  supplied  50/50  ramp  function  that  was  applied  to  the  diode  driver  at  a 
frequency  of  400  Hz  with  a  peak-to-peak  voltage  of  40  mV.  These  settings  ensured 
the  diode  would  continuously  scan  forwards  and  backwards  across  the  rubidium  Di 
absorption  spectrum.  The  diode  temperature  controller  was  altered  to  view  various 
areas  of  the  spectrum  if  desired.  The  signal  generator  trigger  output  was  connected  to 
a  300  MHz  oscilloscope  to  monitor  the  spectrum.  Polarization  beam  splitters  (PBS) 
were  used  to  direct  the  beam  into  four  New  Focus  large-area  infrared  photoreceivers 
(2033),  as  seen  in  Figure  8.  A  half-wave  plate  centered  at  the  diode  wavelength  was 
inserted  before  each  PBS  to  adjust  the  polarization  of  the  beam.  By  adjusting  the 
polarization  with  the  half-wave  plate,  the  amount  of  power  (<  20  mW )  projected 
on  each  detector  could  be  modified  to  avoid  transient  broadening  and  saturating  the 
photoreceivers.  The  half-wave  plates  were  turned  such  that  identical  signals  were 
being  displayed  on  the  oscilloscope  using  the  isolated  beam  with  no  spectrum  being 
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Figure  8.  Experimental  apparatus  used  to  observe  if  rubidium  vapor 
diffused  through  the  capillary  tube  into  the  pump  side  chamber.  The 
red  lines  designate  the  Di  diode  laser  and  the  green  lines  represent  the 
D2  diode  laser  (the  color  of  the  lines  are  not  reflective  on  the  frequency 
of  the  lasers  but  simply  to  differentiate  where  each  laser  was  directed). 
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collected.  A  low  pressure  rubidium  reference  cell  containing  8'Rb  was  used  to  con¬ 
firm  the  laser  scanning  frequency  and  assist  in  finding  the  Di  spectrum.  The  detector 
receiving  the  transmitted  signal  from  the  reference  cell,  Iref,  was  connected  to  the  os¬ 
cilloscope.  A  second  PBS  directed  the  scanning  beam  into  the  fore  pressure  chamber 
containing  both  isotopes  of  rubidium  in  their  natural  abundances.  The  light  trans¬ 
mitted  through  the  fore  pressure  chamber  was  collected  by  a  different  photoreceiver, 
Ii .  Another  PBS  was  introduced  guiding  the  laser  into  the  back  pressure  chamber  in 
order  to  reveal  if  any  rubidium  successfully  diffused  through  the  tube  and  into  the 
pump  side  chamber  using  a  detector,  J2.  A  fourth  photoreceiver,  J0,  captured  the 
initial  scanning  laser  light  and  was  used  to  determine  the  alkali  vapor  number  density 
through  Beer’s  Law  in  equation  (37). 

Two  experiments  were  conducted  using  two  different  capillary  tube  lengths.  The 
first  tube  used  was  20  cm  long  and  the  ends  of  the  tube  were  positioned  approximately 
1  cm  from  the  edges  of  each  chamber,  displayed  in  Figure  9a.  As  shown  in  Figure  9b, 
to  avoid  clipping  the  capillary  tube  as  the  laser  propagated  through  the  chambers, 
the  beam  was  directed  above  the  tube.  The  second  experiment  used  a  12  cm  long 
capillary  tube  and  the  ends  were  located  approximately  in  the  center  of  each  chamber, 
displayed  in  Figure  10a.  For  this  case,  the  beams  were  directed  approximately  0.5  cm 
away  from  the  ends  of  the  tube,  shown  in  Figure  10b.  A  Tiger  diode  laser  tuned  at 
the  rubidium  D2  resonate  frequency  (approximately  780  nm )  with  a  beam  diameter 
of  approximately  0.25  cm  was  aligned  through  the  500  /im  inner-diameter  capillary 
tube  using  two  highly-reflective  (HR)  mirrors.  Because  of  the  small  diameter  of  the 
capillary  tube,  only  about  4%  of  the  beam  was  coupled  into  the  tube.  The  400  mW 
diode  source  acted  as  a  possible  ablation  mechanism  for  any  rubidium  atoms  inter¬ 
acting  with  pre-existing  contaminants  formulating  on  the  inner-walls  of  the  capillary 
tube.  A  higher  power  ablation  source  is  commonly  utilized,  but  the  Tiger  diode  laser 
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Figure  9.  (a)  20  cm  long  capillary  tube  location  with  both  ends  approx¬ 
imately  1  cm  from  the  edges  of  the  chambers,  (b)  Location  of  the  Di 
scanning  laser  beam  in  relation  to  the  capillary  tube. 
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Figure  10.  (a)  12  cm  long  capillary  tube  location  with  both  ends  located 
near  the  center  of  the  chambers,  (b)  Location  of  the  Di  scanning  laser 
beam  in  relation  to  the  capillary  tube. 
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satisfied  the  tuning  and  readily  availability  requirements.  The  two  HR  mirrors  were 
fine-tuned  until  the  majority  of  the  light  propagated  down  the  tube  with  minimal 
transmission  losses. 
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IV.  Results  and  Analysis 


Flow  Characterization  and  Viscosity 

The  viscosity  of  a  gas  confirms  the  type  of  flow  through  the  capillary  tube.  Dif¬ 
fusion  rates  for  each  gas  exploit  important  pressure  differential  limits  that  must  be 
considered.  Remaining  within  the  boundaries  will  ensure  laminar  flow  and  avoid 
jeopardizing  the  gaseous  state  of  the  rubidium.  Once  a  controlled  flow  environment 
is  implemented,  velocity  profile  distributions  can  be  formed.  The  flow  profiles  re¬ 
veal  which  gas  propagates  through  the  capillary  tube  with  the  largest  velocity  while 
remaining  within  the  constraints  of  Poiscuille  flow  theory. 

In  four  separate  experiments,  approximately  600  Torr  of  a  noble  gas  was  added 
to  the  fore  pressure  chamber.  As  the  gas  diffused  through  the  glass  rod,  the  pressure 
decrease/increase  of  both  chambers  was  recorded  by  the  two  pressure  sensors,  as  seen 
in  Figure  11.  The  diffusion  rates,  shown  in  Figure  12,  contain  the  decreasing  pressure 
in  the  fore  chamber  for  each  of  the  four  noble  gases  used  in  the  experiment.  In  Fig¬ 
ure  12,  the  data  collection  began  where  the  fore  pressure  reached  600  Torr  so  they 
would  be  approximately  on  the  same  initial  pressure  scale.  Where  the  cubes  came 
to  equilibrium  varied  due  to  varying  starting  pressures  and  the  rate  at  which  the  gas 
was  added  to  the  fore  pressure  chamber.  Figure  13  displays  the  reciprocal  of  the  fore 
pressure  data,  which  is  used  to  extract  the  viscosity  of  the  gas.  This  is  the  first  step 
in  the  analysis  to  determine  the  viscosities  of  each  gas  as  it  flows  through  the  glass 
rod.  Once  the  viscosities  are  calculated  for  the  volumetric  parameters  of  the  system, 
assumptions  can  be  made  about  flow  velocities  and  the  associated  propagating  profile 
distributions.  These  characteristics  are  important  in  determining  which  gas  would  be 
most  suitable  to  be  used  as  a  potential  carrier  of  rubidium  vapor  atoms  through  the 
capillary  tube.  In  the  initial  seconds  of  the  inverse  pressure  rates,  where  the  slopes  ap- 
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Figure  11.  Pressure  readings  from  the  fore  and  back  pressure  chambers 
as  approximately  600  Torr  of  krypton  was  added  to  the  fore  pressure 
cube  and  diffused  to  the  back  pressure  chamber.  Where  the  cubes  came 
to  equilibrium  varied  due  to  varying  starting  pressures  and  the  rate  at 
which  the  gas  was  added  to  the  fore  pressure  chamber. 


Figure  12.  Fore  pressure  decrease  data  as  the  noble  gases  diffused 
through  the  glass  rod  to  the  back  pressure  chamber.  Each  gas  data 
set  starts  at  approximately  600  Torr  in  order  to  have  the  same  relative 
initial  pressure  scale. 
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Figure  13.  The  inverse  fore  pressure  data  for  each  of  the  used  noble 
gases. 


pear  to  be  nearly  linear,  the  data  was  fit  to  a  basic  y  =  mx  +  b  linear  equation  as  seen 
in  Figure  14.  The  linear  regime  was  selected  where  the  residual  was  nearly  zero  over 
the  longest  period  of  time.  Least-Squares  linear  regression  provided  minimal  resid¬ 
ual  to  the  fit  compared  to  the  attempted  Levenberg-Marquardt  and  Gauss-Newton 
methods.  The  fit  with  residuals  are  shown  in  Figure  15.  The  slopes  were  inserted 
into  equation  (19)  as  the  time  dependant  derivative  of  the  inverse  fore  pressure  and 
the  viscosities  for  each  noble  gas  were  extracted  and  compared  to  previously  recorded 
values  by  NIST,  as  seen  in  Table  2. 

Establishing  the  viscosities  by  fitting  the  linear  portion  of  the  experimental  inverse 
fore  pressure  data  was  an  approach  used  in  previous  literature.  The  linear  fitting 
technique  assumes  p\  3>  p2.  However,  including  the  back  pressure  measurements 
considers  all  the  physical  properties  occurring  on  both  ends  of  the  capillary  tube  and 
provides  more  information  about  the  flow,  equation  (22)  explores  all  of  the  fluid 
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Figure  14.  Inverse  fore  pressure  data  was  fit  to  a  y  =  mx+b  linear  equation 
for  krypton.  The  slope  of  the  fit  equation  was  used  to  determine  the 
viscosity  of  the  flowing  gas  within  the  capillary  tube.  The  data  inside 
the  dashed  box  represents  the  region  that  was  appeared  nearly  linear. 
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Figure  15.  (a)  Linear  fit  to  krypton  experimental  inverse  fore  pressure 
data  using  the  Least-Squares  linear  regression  method,  (b)  Residual 
between  the  raw  data  and  the  best  fit  line. 
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Table  2.  The  viscosities  for  fitting  the  slope,  fitting  the  entire  curve  and 
the  previous  results  recorded  by  NIST  are  shown  in  fiPa  —  Sec. 


Noble  Gas 

Slope  Fit 

Curve  Fit 

NIST  Results 

Helium 

61.310  ±  1.566 

5.420  ±  0.238 

19.838  ±1.984 

Argon 

53.254  ±  1.213 

20.604  ±0.251 

22.543  ±2.254 

Krypton 

73.405  ±2.187 

22.912  ±0.259 

25.111  ±2.511 

Xenon 

82.934  ±0.919 

36.719  ±0.281 

22.945  ±2.295 

dynamics  occurring  in  the  system  for  any  pressure  differential.  The  fact  that  the 
back  pressure  began  to  increase  before  the  gas  inflow  was  closed  off  was  neglected  in 
the  linear  fit  analysis.  However,  the  back  pressure  employs  a  significant  role  when 
fitting  the  entire  inverse  fore  pressure  curve.  The  definition  of  the  total  pressure  in 
the  system  is  sum  of  the  pressures  contained  in  both  chambers.  When  the  gas  began 
entering  the  fore  pressure  chamber,  the  diffusion  process  had  already  commenced 
because  there  was  no  shut-off  valve  between  the  two  chambers.  To  accommodate  this 
situation,  the  initial  back  pressure  was  subtracted  from  all  of  the  fore  pressure  values 
confirming  equation  (21).  From  the  solution  to  equation  (24),  the  Least-Squares 
method  was  utilized  to  fit  the  entire  experimental  inverse  fore  pressure  curves  with 
two  fit  parameters.  A  demonstrated  fit  along  with  residuals  between  the  data  and  the 
fit  is  shown  in  Figure  16.  The  fit  parameter,  T,  was  then  substituted  into  equation 
(25)  and  the  viscosities  for  each  gas  were  determined.  The  viscosities  for  the  linear  fit 
and  the  ODE  curve  fit  were  converted  to  the  commonly  published  units  of  fiPa  —  Sec, 
as  seen  in  Table  2.  By  fitting  the  entire  ODE  characterizing  all  the  fluid  dynamics 
of  the  system  rather  than  exclusively  the  linear  segment  of  the  inverse  fore  pressure 
yields  viscosity  measurements  closer  to  the  previously  published  values. 

In  order  to  understand  the  type  of  flow,  laminar  or  turbulent,  occurring  within 
the  glass  rod,  the  Reynolds  Number  must  be  determined.  Back  pressures  for  each 
experiment  were  subtracted  from  the  fore  pressures  to  find  the  actual  pressure  gradient 
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Figure  16.  (a)  ODE  curve  fit  to  krypton  experimental  inverse  fore  pres¬ 
sure  data  using  the  Least-Squares  linear  regression  method,  (b)  Residual 
between  the  raw  data  and  the  best  fit  curve. 


at  each  point  in  time.  This  determines  the  maximum  flow  velocity  in  equation  (30), 
which  will  assist  in  choosing  a  carrier  gas.  Using  the  density  of  each  gas  at  room 
temperature  and  resolved  viscosities,  the  Reynolds  Numbers  can  be  calculated  with 
equation  (27).  Figure  17  and  Figure  18  display  the  Reynolds  Number  for  each  gas 
at  every  pressure  differential  for  the  viscosities  extracted  from  the  linear  fit  and  the 
ODE  fit,  respectively. 

If  the  Reynolds  Number  is  less  than  approximately  30,  then  the  flow  through  the 
tube  can  be  assumed  to  follow  Poiseuillc  flow  theory.  When  the  Reynolds  Number  is 
greater  than  30,  this  simply  means  the  gas  does  not  enter  the  tube  with  a  parabolic 
velocity  profile  described  in  equation  (7).  If  the  flow  is  turbulent,  determining  the  flow 
distribution  profiles  and  velocities  becomes  a  non-trivial  procedure.  At  some  length 
down  the  tube,  the  flow  distribution  may  develop  into  a  Poiseuille  flow  model  and 
the  velocity  profiles  can  be  extracted.  The  length  down  the  tube  from  the  gas  entry 
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Figure  17.  Left  Axis:  Reynolds  Number  at  each  pressure  differential  for 
each  of  the  noble  gases  using  the  linear  fit  viscosity  value.  Right  Axis: 
Length  down  the  12  cm  tube  before  Poiseuille  flow  theory  can  be  applied 
and  the  parabolic  velocity  profile  is  fully  developed. 


Figure  18.  Left  Axis:  Reynolds  Number  at  each  pressure  differential  for 
each  of  the  noble  gases  using  the  ODE  curve  fit  viscosity  value.  Right 
Axis:  Length  down  the  12  cm  tube  before  Poiseuille  flow  theory  can  be 
applied  and  the  parabolic  velocity  profile  is  fully  developed. 
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point  that  Poiseuille  theory  applies  is  found  by  using  equation  (29).  From  the  results 
shown  in  Figure  17  following  the  linear  fit  method,  it  can  be  seen  that  Poiseuille 
flow  theory  is  only  present  within  the  tube  for  small  pressure  differentials  (A p  <  10 
Torr ).  The  flow  is  then  considered  to  be  turbulent  for  0.25  to  almost  3  cm  down 
the  capillary  for  larger  pressure  differentials  (with  the  exception  of  helium,  which 
remains  in  the  Poiseuille  flow  regime  up  to  a  146  Torr  pressure  differential).  The 
length  of  the  capillary  was  measured  to  be  12  cm  long,  revealing  25%  of  the  tube 
is  under  unpredictable  and  chaotic  turbulent  flow.  The  viscosities  determined  by 
fitting  the  ODE  curve  produced  Reynolds  Numbers  that  were  approximately  a  factor 
of  two  smaller  than  the  ones  calculated  through  the  linear  fit  method,  as  seen  in 
Figure  18.  These  Reynolds  Numbers  convey  a  larger  pressure  differential  margin  can 
be  explored  before  turbulent  flow  is  experienced.  The  accelerated  routine  of  solely 
fitting  the  initial  inverse  pressure  data  produces  larger  viscosity  values  and  in  turn 
outputs  larger  Reynolds  Numbers.  The  technique  explored  in  this  analysis  achieved 
more  precise  viscosity  measurements  yielding  lower  Reynolds  Numbers  describing  the 
flow  through  the  capillary  tube.  Manifesting  the  differential  equations  describing  the 
known  dynamics  of  the  system  and  formulating  solutions  that  mapped  experimental 
data  proved  to  be  a  justifiable  avenue  and  modeling  the  drastically  chaotic  turbulent 
behavior  was  avoided. 

When  considering  techniques  on  transferring  number  densities  of  alkali  vapor 
through  the  tube  used  in  this  experiment,  one  such  technique  that  can  be  consid¬ 
ered  is  adding  one  of  the  noble  gasses  and  allowing  diffusion  to  transfer  the  alkali. 
Maintaining  a  pressure  differential  between  the  two  chambers  less  than  150  Torr 
permits  an  accurate  description  of  the  flow  velocity  (without  any  alkali)  through  the 
tube.  The  velocity  flow  distributions  under  Poiseuille  flow  theory  coincide  with  equa¬ 
tion  (7)  where  the  profile  of  the  gas  propagating  through  the  tube  can  be  described  by 
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Figure  19.  Velocity  profile  distribution  for  helium  as  a  function  of  dis¬ 
tance  from  the  center  of  the  glass  rod,  r,  and  for  the  pressure  differentials 
that  follow  Poiseuille  flow  theory  determined  by  the  Reynolds  Number 
being  less  than  30. 


a  parabola.  The  flow  velocity  gradient  is  determined  from  equation  (7)  as  a  function 
of  r,  the  distance  from  the  center  of  the  tube.  Figure  19  through  Figure  22  maps 
out  the  velocity  profiles  for  each  gas  as  a  function  of  distance  from  the  center  of  the 
tube  only  when  the  flow  through  the  tube  is  described  by  100%  Poiseuille  flow  theory 
determined  from  fitting  the  ODE  as  displayed  in  Figure  18. 

Transferring  alkali  through  small  inner-diameter  tubes  has  been  considered  to  be 
a  time  consuming  procedure  because  previous  techniques  use  basic  diffusion  (with  no 
pumping  or  carrier  gases  used)  and  work  on  a  time  scale  of  weeks  or  even  months 
to  provide  sufficient  results.  By  applying  a  noble  gas  to  assist  in  the  propagation 
process,  the  window  of  time  can  be  decreased  dramatically.  From  the  flow  velocity 
results,  helium  propagates  with  the  largest  velocity  through  the  tube  with  a  flow  that 
can  be  fully  understood  for  a  restricted  pressure  differential  regime. 
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Figure  20.  Velocity  profile  distribution  for  argon  as  a  function  of  distance 
from  the  center  of  the  glass  rod,  r,  and  for  the  pressure  differentials  that 
follow  Poiseuille  flow  theory  determined  by  the  Reynolds  Number  being 
less  than  30. 
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Figure  21.  Velocity  profile  distribution  for  krypton  as  a  function  of  dis¬ 
tance  from  the  center  of  the  glass  rod,  r,  and  for  the  pressure  differentials 
that  follow  Poiseuille  flow  theory  determined  by  the  Reynolds  Number 
being  less  than  30. 
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Figure  22.  Velocity  profile  distribution  for  xenon  as  a  function  of  distance 
from  the  center  of  the  glass  rod,  r,  and  for  the  pressure  differentials  that 
follow  Poiseuille  flow  theory  determined  by  the  Reynolds  Number  being 
less  than  30. 


Rubidium  Number  Densities 

First  attempts  to  transfer  rubidium  vapor  through  the  500  /irn  capillary  tube 
consisted  of  pumping  on  the  back  pressure  chamber  with  the  turbo  pump.  The  system 
was  heated  to  150  C  to  generate  a  rubidium  vapor  in  the  fore  pressure  cube.  A  carrier 
gas  was  not  introduced  for  this  procedure  with  the  anticipation  that  rubidium  atoms 
would  propagate  through  the  tube  with  60  nTorr  in  the  aft  chamber.  The  capillary 
tube  used  was  20  cm  in  length  and  the  ends  of  the  tubes  were  approximately  1  cm 
from  the  edges  of  the  cubes,  shown  previously  in  Figure  9a.  Scanning  diode  laser 
beams  were  directed  slightly  above  the  capillary  tube  searching  for  the  rubidium  Di 
absorption  spectrum,  displayed  in  Figure  9b.  A  significant  spectrum  was  extracted 
on  the  fore  pressure  chamber  when  the  system  was  heated,  but  there  was  no  trace 
of  a  spectrum  on  the  back  pressure  chamber.  Recognizing  the  time  constraints  of 
the  project,  an  alternate  technique  was  deployed  after  five  days  of  pumping  with  no 
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apparent  visible  spectrum. 

Importing  a  carrier  gas  was  the  next  technique  examined.  Heat  tape  was  wrapped 
around  the  main  bus  importing  the  noble  gases  to  avoid  changing  the  temperature 
within  the  chamber  system.  Approximately  40  to  80  Torr  of  krypton  was  added  to 
the  fore  pressure  chamber  with  no  pumping  mechanism  on  the  back  pressure  chamber. 
After  the  gas  diffused  through  the  capillary  into  the  back  pressure  cube,  the  gas  was 
slowly  released  to  the  pumping  system  and  the  process  was  repeated.  On  the  fifth 
day  of  practicing  this  approach,  the  fore  pressure  chamber  held  a  constant  pressure 
of  40  Torr  while  the  back  side  cube  maintained  a  pressure  of  a  few  mTorr.  This 
proposed  an  apparent  blockage  within  the  capillary  and  discontinuing  the  efforts  of 
diffusing  rubidium  vapor  with  a  carrier  gas.  The  entire  system  was  disassembled  to  be 
disinfected  and  rebuilt.  During  the  procedure  of  dismantling  the  apparatus,  a  visible 
magnitude  of  rubidium  accumulation  was  discovered  on  the  calcium  fluoride  window 
nearest  to  the  exit  of  the  tube  in  the  back  pressure  chamber,  as  seen  in  Figure  23. 
The  deposit  formed  a  1  cm  diameter  circle  on  the  window  and  developed  a  raised 
conical  formation  in  the  center  with  a  FWHM  of  approximately  0.05  cm. 

Immediately  after  removing  the  capillary  from  the  chamber  the  mass  was  recorded 
to  be  13.899  grams.  Directing  a  white  LED  source  axially  down  the  length  of  the 
tube  exposed  rubidium  formations  from  the  scattering  light  (not  from  fluorescence), 
displayed  in  Figure  24  and  Figure  25.  Ethanol  was  injected  into  the  tube  cleansing 
it  of  nearly  all  inhabiting  contaminants  and  was  then  dried  by  heat  enforcement.  The 
mass  of  the  evacuated  tube  was  measured  again  at  13.893  grams,  revealing  approx¬ 
imately  6  milligrams  of  substance  within  the  tube,  which  makes  up  9.92%  of  total 
volume  of  the  tube.  This  evidence  proposes  that  as  the  rubidium/krypton  mixture 
propagated  through  the  tube,  atoms  were  adhering  to  the  inner  walls  of  the  capil¬ 
lary.  If  the  rubidium  volume  to  tube  volume  ratio  increases  over  time,  yielding  more 
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Figure  23.  Visible  rubidium  accumulation  (with  an  outer  diameter  of  1 
cm  with  a  raised  cone  in  the  center  having  a  FWHM  of  approximately  0.05 
cm)  on  the  calcium  fluoride  window  closest  to  the  exit  of  the  capillary  on 
the  pump  side  chamber.  The  end  of  the  tube  was  positioned  about  1  cm 
away  from  the  window  and  is  shown  in  Figure  9.  This  validates  rubidium 
successfully  transferred  through  the  20  cm  long,  500  /im  inner-diameter 
capillary  tube. 


Figure  24.  A  white  LED  source  was  directed  down  the  length  of  the  cap¬ 
illary  tube  exposing  rubidium  deposits  within  the  tube  through  scatter 
light  effects  (not  fluorescence). 
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Figure  25.  Color  inverted  image  of  Figure  24  displaying  rubidium  de¬ 
posits  within  the  capillary  through  scatter  light  effects  (not  fluorescence) 
from  a  white  LED  source. 

atomic  layers  condensing  to  the  walls  of  the  tube,  the  flow  rates  will  decrease.  This  is 
because  the  volumetric  flow  rate  is  proportional  to  i?4,  shown  in  equation  (10),  which 
will  make  the  task  of  transferring  rubidium  near  impossible.  The  agents  responsible 
for  the  agglomeration  were  postulated  to  be  derivatives  of  oxidized  substances  resid¬ 
ing  within  the  tube  after  the  manufacture  process.  Decontamination  procedures  were 
exerted  prior  to  the  experiments  by  heating  and  pumping  down  the  system  for  several 
days.  Further  cleansing  methods  should  be  formulated  in  future  work  to  prevent  such 
circumstances. 

New  approaches  were  developed  striving  for  a  successful  rubidium  absorption  spec¬ 
trum.  The  length  of  the  capillary  tube  was  decreased  to  12  cm  and  the  ends  of  the 
tube  were  localized  to  the  center  of  each  chamber  rather  than  close  proximity  to  the 
edges  of  the  cubes,  shown  in  Figure  10a.  This  would  decrease  the  diffusion  time.  Also, 
with  the  ends  of  the  tube  in  the  center  of  the  pressure  chambers,  the  scanning  laser 
sources  could  be  directed  approximately  0.5  cm  from  the  entrance  and  exit  holes  of  the 
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capillary,  shown  in  Figure  10b.  Probing  the  endpoints  of  the  tube  allow  a  spectrum 
to  be  acquired  from  any  alkali  vapor  flowing  through  the  tube.  From  the  previous  two 
techniques,  exclusively  pumping  with  no  carrier  gas  and  basic  diffusion,  the  remedy 
for  successfully  transferring  rubidium  through  the  tube  was  realized.  From  the  fluid 
dynamical  analysis,  helium  proved  to  implement  the  highest  flow  velocity  with  a  vast 
permissible  pressure  differential  margin  while  being  constrained  by  the  Poiseuille  flow 
limits.  Remaining  within  the  Poiseuille  flow  regime  ensured  no  chaotic  flow  was  oc¬ 
curring  that  could  produce  phenomena  like  shock  waves  that  might  solidify  the  alkali 
in  the  capillary.  With  the  vacuum  pump  slightly  opened,  a  needle  valve  controlling 
the  warm  helium  source  was  turned  such  that  a  slow  steady  constant  flow  of  gas  was 
being  introduced  to  the  system.  A  constant  pressure  gradient  of  approximately  80 
Torr  was  maintained  throughout  this  procedure.  Continuously  flowing  the  rubidi- 
um/helium  combination  through  the  system  furnished  fruitful  results.  An  absorption 
spectrum  representing  the  Dj  atomic  hyperhne  transitions  for  both  isotopes  of  rubid¬ 
ium  were  revealed,  shown  in  Figure  26.  Changing  the  pressure  gradient  across  the 
capillary  increased  the  flow  of  alkali  vapor  resulting  in  larger  number  densities  and 
the  peaks  became  even  more  evident  as  seen  in  Figure  27  and  Figure  28  for  pressure 
gradients  of  approximately  120  and  150  Torr,  respectively. 

In  order  to  determine  the  number  density  of  the  alkali  vapor  that  transferred 
into  the  back  pressure  chamber,  a  summation  of  Voigt  profiles  was  developed  to  fit 
the  four  hyperhne  transitions  for  each  of  the  isotopes  using  equation  (46).  If  the 
laser  probe  has  an  intensity  greater  than  the  saturation  intensity  ( I/Isat  >  1),  the 
incident  intensity  is  considered  a  pump  source  rather  than  a  probe.  Because  of  this 
condition,  atoms  excited  to  an  upper  hyperhne  component  may  relax  to  a  different 
hyperhne  ground  state  energy  level  than  which  it  initially  started  due  to  spontaneous 
or  stimulated  emission.  For  example,  if  the  F”  =  2  — y  F'  =  3  transition  in  85Rb  is 
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Figure  26.  Atomic  rubidium  absorption  peaks  being  collected  from  the 
back  pressure  chamber.  Results  accomplished  by  maintaining  a  constant 
pressure  gradient  of  80  Torr  across  the  capillary  tube  with  a  constant 
input  flow  of  helium  in  the  fore  pressure  chamber  and  the  vacuum  pump 
slightly  opened  in  the  back  pressure  chamber. 
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Figure  27.  Atomic  absorption  peaks  being  collected  from  the  back  pres¬ 
sure  chamber.  A  pressure  gradient  of  approximately  120  Torr  was  estab¬ 
lished  to  increase  the  number  density  flowing  through  the  capillary. 


Figure  28.  Atomic  absorption  peaks  being  collected  from  the  back  pres¬ 
sure  chamber.  A  pressure  gradient  of  approximately  150  Torr  was  estab¬ 
lished  to  increase  the  number  density  flowing  through  the  capillary. 
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Figure  29.  Fit  to  the  multi-isotopic  rubidium  absorbance  spectrum  that 
was  transferred  through  the  capillary  tube  and  into  the  pump  chamber 
along  with  residuals  displaying  the  difference  between  the  fit  (red)  and 
the  raw  data  (blue)  with  A P  =  120 Torr. 


optically  pumped,  the  atom  may  relax  to  the  lower  F"  =  3  ground  state  rather  than 
the  F"  =  2  state.  By  depleting  the  hyperfine  ground  state  being  investigated  directly 
effects  the  predicted  Boltzmann  distribution  in  equation  (44).  This  will  change  the 
relative  amplitude  of  the  interrogated  transition  and  the  practice  of  fitting  the  Voigt 
profiles  will  become  difficult.  From  the  spectrum  fit  in  Figure  29,  the  residual  is 
calculated  by  subtracting  the  fit  from  the  data.  The  fit  amplitudes  are  larger  than  the 
raw  data  amplitudes  only  for  the  87Rb  hyperfine  structure,  but  are  smaller  than  the 
data  for  the  85 Rb  spectrum  (see  fit  and  residual  in  Figure  29.  It  can  be  inferred  that 
the  isotopic  ratio  has  changed  from  the  naturally  occurring  abundances  previously 
stated.  The  8'Rb  atom  has  a  larger  mass  than  the  85 Rb  atom.  Because  of  the  mass 
difference  the  heavier  atoms  have  different  flow  characteristics  than  the  lighter  atoms. 
To  estimate  number  densities  with  more  precision,  further  analysis  must  be  conducted 
to  investigate  the  isotopic  ratio  shift.  For  the  case  of  this  experiment,  the  number 
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densities  will  only  be  slightly  effected  and  the  focus  was  to  approximate  a  number 
density  for  fiber  laser  applications. 

The  absorbance  is  dependant  on  the  path  length  of  the  test  cell.  Since  this  ex¬ 
periment  demonstrates  a  constant  flow  of  alkali  rather  than  a  cell  under  equilibrium 
conditions,  the  path  length  depends  on  the  distance  away  from  the  end  of  the  tube 
the  probe  is  investigating.  The  flow  of  rubidium  exiting  the  capillary  is  modeled  by  a 
constant  diverging  cone.  The  FWHM  of  the  deposit  on  the  window  is  approximately 
0.05  cm  in  diameter.  Determining  an  interaction  path  length  from  this  criteria  was 
determined  for  three  locations  by  mapping  a  constant  diverging  cone  from  the  exit 
of  the  tube  to  the  deposit.  The  first  location  dissected  the  cone  immediately  after 
the  rubidium  exited  the  tube  where  the  path  length  equals  the  0.05  cm  diameter  of 
the  tube.  The  second  location  of  0.15  cm  was  determined  from  the  FWHM  of  the 
rubidium  deposit  on  the  window  of  the  cell.  The  probe  was  aimed  1.0  cm  away  from 
the  end  of  the  tube  which  is  the  same  distance  the  tube  was  when  the  rubidium 
accumulated  on  the  window.  The  third  case  used  the  width  of  the  chamber,  7.0  cm, 
investigating  the  maximum  distance  away  from  the  tube. 

The  Gaussian  and  Lorentzian  lineshape  widths  were  permitted  to  vary  in  the  fit 
because  the  I/Isat  condition  altered  the  population  of  the  hyperhne  ground  states 
which  directly  effects  the  Boltzmann  distribution.  From  Figure  29,  the  constant 
K 0  represents  the  total  amplitude  of  the  profile,  K 1  represents  the  frequency  shift, 
K 2  represents  the  Lorentzian  width,  K 3  represents  the  Gaussian  width,  and  the 
constants  K 4  through  K 6  represent  a  quadratic  baseline  on  either  side  of  the  profile 
distribution  for  a  better  fit  ((K 4)z/2  +  (K5)u  +  (K 6)).  The  absorbance  amplitude  is 
determined  by  a  normalization  factor  of  s/k I \/hi2  which  is  multiplied  by  the  profile 
amplitude  fit  parameter  (K 0).  The  number  density  is  calculated  by  dividing  the 
absorbance  amplitude  by  (/Aq)/(87tt2i).  The  calculated  number  densities  associated 
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with  path  length  are  shown  in  Table  3.  A  rubidium  vapor  pressure  can  be  determined 


Table  3.  Calculated  rubidium  number  densities  for  the  absorbance  data 


fit  in  Figure  29  with  A P  =  YlQTorr  and  estimated  path  lengths. 


Path  Length  (cm) 

Number  Density  (cm  3) 

Error  (cm  3) 

0.05 

1.384  657  x  1012 

2.796  86  x  1010 

0.15 

4.615  525  x  1011 

3.342  76  x  109 

7.00 

9.890  410  x  109 

3.135  48  x  107 

using  the  calculated  number  densities  for  the  three  interaction  path  lengths  previously 
discussed.  The  back  chamber  maintained  a  total  pressure  of  20  Torr.  A  mixing  ratio 
was  determined  by  dividing  the  rubidium  vapor  pressure  by  the  total  pressure  in 
the  chamber.  The  mixing  ratios  are  shown  in  Table  4.  A  rubidium  vapor  pressure 


Table  4.  Mixing  ratios  between  the  rubidium  and  helium  combinations 
are  given  for  the  three  interaction  path  lengths  which  determined  the 
rubidium  number  densities  that  propagated  through  the  capillary  tube. 


Path  Length  (cm) 

Vapor  Pressure  (Torr) 

Mixing  Ratio 

0.05 

5.064  x  10"5 

2.532  ppm 

0.15 

1.688  x  10”5 

0.844  ppm 

7.00 

3.617  x  10“7 

18.085  ppb 

was  estimated  to  be  5.664  x  10~5  Torr  using  the  temperature  of  the  chamber,  which 
was  353.15  K .  The  back  pressure  chamber  maintained  a  total  pressure  of  140  Torr, 
providing  a  gas  mixing  ratio  of  405  ppb.  The  smaller  two  interaction  path  lengths 
from  Table  4  do  not  accurately  quantify  the  mixing  ratio  throughout  the  back  pressure 
chamber.  However,  the  length  of  7.00  cm  is  equivalent  to  the  width  of  the  chamber 
and  can  accurately  describe  the  vapor  pressure.  When  comparing  the  mixing  ratio 
of  the  fore  and  back  pressure  chambers,  it  can  be  seen  that  the  ratio  on  the  fore 
pressure  side  is  approximately  20  times  larger  than  the  gas  ratio  in  the  back  pressure 
chamber.  This  concludes  that  a  sufficient  amount  of  rubidium  is  propagating  through 
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the  capillary  tube  to  the  back  pressure  chamber  without  condensing  to  the  inner  walls 
to  the  point  where  the  radius  of  the  tube  is  effecting  the  flow  of  the  rubidium. 
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V.  Conclusions 


Fluid  Dynamic  Analysis 

This  work  provided  a  thorough  flow  characterization  of  four  noble  gases  to  de¬ 
termine  the  gas  viscosity  and  velocity  profiles  distributions.  Previous  techniques  for 
extracting  flow  viscosities  were  explored  and  compared  to  an  alternate  approach. 
Considering  the  entire  ordinary  differential  equation  that  contains  the  temporal  back 
pressure  measurements  is  critical  in  modeling  the  physical  properties  of  the  system. 
Once  viscosity  values  were  compared  to  recorded  NIST  values,  Poiseuille  flow  theory 
was  investigated.  Following  the  constraints  of  Poiseuille  flow  allowed  for  a  com¬ 
plete  understanding  of  the  velocity  profile  distribution  throughout  the  capillary  tube. 
Establishing  the  type  of  flow  and  velocities  each  gas  experienced  through  the  experi¬ 
ments  aided  in  the  selection  process  of  which  atoms  could  be  used  as  carriers  for  the 
rubidium  vapor.  Avoiding  chaotic  turbulent  flow  was  also  valuable  information  as¬ 
certained  from  the  flow  experiments  which  could  have  lead  to  shock  wave  phenomena 
that  could  have  dismantled  the  propagation  flow.  The  Reynolds  Number  calcula¬ 
tions  outlined  a  pressure  differential  margin  that  would  lead  to  a  quicker  transferring 
routine  without  altering  the  system  parameters  and  apparatus. 

Transferring  Rubidium  Vapor  Through  a  Capillary  Tube 

This  work  examined  three  techniques  on  how  to  transfer  atomic  rubidium  vapor 
through  a  small  inner-diameter  capillary  tube.  Pumping  on  the  system  with  a  turbo 
pump  proved  to  be  a  preliminary  approach  with  no  direct  achievement,  but  portrayed 
a  major  role  in  depicting  a  successful  route.  Basic  diffusion  principals  were  followed 
by  adding  various  pressures  of  krypton  to  the  fore  pressure  chamber  and  allowing  the 
system  to  come  to  equilibrium.  An  absorption  spectrum  was  not  extracted,  but,  a 
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visible  conglomerate  was  evident  on  the  window  closest  to  the  exit  of  the  capillary.  6 
milligrams  of  rubidium  were  removed  from  within  the  tube  revealing  evidence  that 
a  fraction  of  the  propagating  rubidium  was  adhering  to  the  walls  of  the  capillary. 
This  could  have  developed  through  previous  contamination  of  oxidation  molecules 
inhabiting  the  tube  from  the  time  of  manufacture  to  the  time  of  implementation 
in  the  transferring  apparatus.  A  combination  of  the  two  previous  approaches  was 
attempted  and  a  Di  rubidium  absorption  spectrum  was  obtained.  From  the  measured 
spectrum,  number  densities  of  9.89  x  1010  cm-3  and  5.5  x  1010  cm"3  were  calculated  for 
a  path  length  of  7  cm.  DPALs  typically  operate  with  number  densities  on  the  order 
of  1013  cm~3,  which  is  much  larger  than  demonstrated  here.  These  are  preliminary 
techniques  for  advancing  alkali  vapor  filled  hollow-core  fiber  laser  systems  and  new 
approaches  are  currently  being  explored. 

Future  Work 

Preparing  the  capillary  tube  for  the  experiment  is  an  essential  step  when  using 
atomic  alkali.  Alkali  is  highly  reactive  with  any  form  of  oxidation  yielding  emphases 
on  decontamination  procedures  to  evacuate  the  inner  walls  of  the  capillary.  Con¬ 
tamination  in  an  alkali  cell  spreads  quickly  bringing  and  abrupt  halt  to  any  existing 
experiments.  Exploring  what  fraction  of  alkali  naturally  clings  to  the  inner  walls  of 
the  capillary  would  prove  to  be  beneficial  to  ensure  transferring  techniques  can  be 
examined  over  long  periods  of  time  without  the  diameter  of  the  tube  closing  and  re¬ 
sisting  vapor  flow.  Adapter  fittings  connecting  various  parts  of  the  entire  apparatus 
should  be  selected  on  the  basis  that  they  can  be  heated  to  high  temperatures  and 
cooled  multiple  times  without  generating  a  leak  within  the  system.  Varying  the  tem¬ 
perature  could  be  effective  if  the  option  of  solidifying  the  rubidium  before  evacuating 
the  cell  is  preferred. 
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Measuring  the  exiting  flow  characteristics  at  the  end  of  the  tube  would  prove  to 
be  beneficial  in  extracting  number  densities  within  the  chambers.  Using  an  infrared 
camera  to  view  the  fluorescence  given  off  by  the  resonating  rubidium  vapor  atoms 
would  provide  a  more  accurate  path  length  and  more  precisely  predict  the  recorded 
number  densities.  Once  a  spectrum  is  extracted  from  the  back  pressure  chamber  and 
rubidium  vapor  is  apparent  on  both  ends  of  the  capillary,  alkali  could  be  assumed 
to  be  within  the  tube.  The  number  densities  inside  the  tube  are  even  more  critical 
to  the  fiber  laser  community,  especially  for  a  diode  pumped  alkali  fiber  amplifier 
application.  The  capillary  tube  was  solely  employed  as  a  surrogate  to  generate  number 
densities  within  a  small  inner-diameter  tube.  Scaling  down  the  inner-diameter  should 
be  examined  providing  thresholds  on  how  small  the  diameter  can  be  before  alkali  flow 
is  restricted  by  the  tube  and  what  maximum  number  densities  of  various  gases  can 
be  issued  inside  the  tube  for  useful  fiber  laser  applications. 

A  more  precise  measurement  of  the  number  density  could  be  determined  if  the 
laser  probe  is  investigating  hyperfine  transitions  at  an  intensity  that  is  less  than  the 
saturation  intensity.  By  neglecting  this  criteria,  the  hyperfine  ground  state  at  line 
center  is  being  depleted  and  the  Maxwell  Boltzmann  distribution  does  not  correctly 
model  the  relative  amplitudes  of  the  transitions.  The  flow  of  alkali  vapor  through 
a  small  inner  diameter  capillary  tube  should  be  accurately  modeled  for  future  work. 
Determining  if  the  mass  difference  between  isotopes  plays  an  important  role  in  the 
flow  velocity  could  assist  in  understanding  why  the  isotopic  ratios  were  shifted  in  the 
rubidium  spectrum. 
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Appendix  A.  Fitting  Inverse  Pressure  Curves 


The  inverse  pressure  curves  were  fit  to  two  different  equations.  The  linear  equation 
assumes  p\  3>  P2  and  fits  the  linear  regime  of  the  inverse  pressure.  The  ODE  fit 
includes  all  of  the  physical  dynamics  of  the  system  and  does  not  neglect  the  pressure 
on  the  back  chamber.  Figure  30  through  Figure  33  show  the  linear  regime  fits  with  the 
residuals.  Figure  34  through  Figure  37  show  the  fit  to  the  ODE  with  the  residuals. 
Table  5  displays  the  fit  parameters  for  the  linear  regime,  specifically  the  slope  of 
the  line  and  the  y-intercept  that  best  fits  the  data  using  the  Least-Squares  linear 
regression  method.  Table  6  displays  the  fit  constants  that  are  directly  related  to  the 
viscosity  as  well  as  the  estimated  total  pressure  constant. 


Table  5.  Linear  fit  parameters  determined  from  Least- Squares  linear 
regression  method. 


Noble  Gas 

Slope  (( Torr  —  sec )  4) 

y-intercept  (Torr  4) 

Helium 

8.723  x  1(T5 

0.0017 

Argon 

1.008  x  AT4 

0.0017 

Krypton 

7.310  x  AT5 

0.0017 

Xenon 

6.470  x  10“5 

0.0016 

Table  6.  ODE  fit  parameters  determined  from  Least-Squares  linear  re¬ 
gression  method. 


Noble  Gas 

Total  Pressure  (po) 

Constants  (T) 

Helium 

216.344 

0.428300 

Argon 

445.591 

0.232063 

Krypton 

408.249 

0.191210 

Xenon 

479.000 

0.140000 
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Figure  30.  (a)  Linear  curve  fit  to  helium  experimental  inverse  fore  pres¬ 
sure  data  using  Least-Squares  linear  regression  method,  (b)  Residual 
between  the  raw  data  and  the  best  fit  curve. 
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Figure  31.  (a)  Linear  curve  fit  to  argon  experimental  inverse  fore  pres¬ 
sure  data  using  Least-Squares  linear  regression  method,  (b)  Residual 
between  the  raw  data  and  the  best  fit  curve. 
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Figure  32.  (a)  Linear  curve  fit  to  krypton  experimental  inverse  fore 

pressure  data  using  Least-Squares  linear  regression  method,  (b)  Residual 
between  the  raw  data  and  the  best  fit  curve. 


Figure  33.  (a)  Linear  curve  fit  to  xenon  experimental  inverse  fore  pres¬ 
sure  data  using  Least-Squares  linear  regression  method,  (b)  Residual 
between  the  raw  data  and  the  best  fit  curve. 
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Figure  34.  (a)  ODE  curve  fit  to  helium  experimental  inverse  fore  pressure 
data  using  Least-Squares  linear  regression  method,  (b)  Residual  between 
the  raw  data  and  the  best  fit  curve. 
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Figure  35.  (a)  ODE  curve  fit  to  argon  experimental  inverse  fore  pressure 
data  using  Least-Squares  linear  regression  method,  (b)  Residual  between 
the  raw  data  and  the  best  fit  curve. 
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Figure  36.  (a)  ODE  curve  fit  to  krypton  experimental  inverse  fore  pres¬ 
sure  data  using  Least-Squares  linear  regression  method,  (b)  Residual 
between  the  raw  data  and  the  best  fit  curve. 


x  10‘3  (a) 


Figure  37.  (a)  ODE  curve  fit  to  xenon  experimental  inverse  fore  pressure 
data  using  Least-Squares  linear  regression  method,  (b)  Residual  between 
the  raw  data  and  the  best  fit  curve. 
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